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ABSTRACT 

Pseudomonas aeruginosa relies on cell motility and 
ability to form biofilms to establish infections; 
however, the mechanism of regulation remains 
obscure. Here we report that BswR, a xenobiotic 
response element-type transcriptional regulator, 
plays a critical role in regulation of bacterial 
motility and biofilm formation in P. aeruginosa. 
Transcriptomic and biochemical analyses showed 
that BswR counteracts the repressor activity of 
MvaT, controls the transcription of small RNA 
rsmZ and regulates the biogenesis of bacterial 
flagella. The crystal structure of BswR was deter- 
mined at 2.3 A resolution; the monomer comprises 
a DNA-binding domain with a helix-turn-helix motif 
in the N terminus and two helices (a6 and a7) with a 
V-shaped arrangement in the C-terminus. In 
addition to the contacts between the parallel 
helices a5 of two monomers, the two helical exten- 
sions (a6 and a7) intertwine together to form a 
homodimer, which is the biological function unit. 
Based on the result of DNase I protection assay 
together with structural analysis of BswR 
homodimer, we proposed a BswR-DNA model, 
which suggests a molecular mechanism with 
which BswR could interact with DNA. Taken 
together, our results unveiled a novel regulatory 
mechanism, in which BswR controls the motility 
and biofilm formation of P. aeruginosa by 
modulating the transcription of small RNA rsmZ. 



INTRODUCTION 

Pseudomonas aeruginosa is an important opportunistic 
human pathogen that causes acute or chronic infections 
in immunocompromised patients (1). The infectivity of 
P. aeruginosa is associated with its motility and biofilm 
formation (2), which is mutually exclusively switching 
between motile and sessile states on environmental 
conditions. Switching between motility and sessility aids 
the pathogen to survive harsh environmental conditions 
by increasing the efficiency of nutrient acquisition, 
escaping from toxic substances and accessing to favorable 
colonization sites (3,4). Mutants lacking motility and 
biofilm formation showed attenuated virulence in a 
burned mouse model (5), reduced invasion in corneal 
epithelial cells (6) and decreased adhesion to human 
respiratory mucin (7). 

Swarming is one of the types of bacterial motilities 
through which bacterial cells move around to aid 
systemic infection and biofilm formation. It is a complex 
adaptation process in response to various environmental 
cues (8). In P. aeruginosa, the swarming motility is 
intimately associated with flagella assembly and type IV 
pilus (T4P) biogenesis (9). The assembly of flagella 
involves >50 genes that are coordinated by an array of 
regulatory proteins including the master regulator FleQ, 
the sigma factor RpoN, transcription factors FleR and 
FliA (10). Pilus biogenesis, an equally complicated 
process, also involves >40 genes and multiple signal 
transduction pathways (11). These signaling pathways 
include the two-component regulatory systems PilR/PilS 
and AlgR/FimS (12), the global regulator Crc (13), the 
virulence factor regulator Vfr (14) and the chemosensory 
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systems encoded by the pilGHIJK-chpABC gene clusters 
(15,16). 

In addition to these regulatory proteins, the small RNA 
regulatory system RsmA/RsmZ is also known to play 
critical roles in the regulation of bacterial motility and 
biofilm formation (17,18). RsmA, an RNA-binding 
protein, regulates swarming by positively modulating the 
biosynthesis of flagella and T4P (19,20). RsmZ, a 
noncoding regulatory small RNA, is an antagonist of 
the functional RsmA. Overexpression of RsmZ inactivates 
RsmA and abolishes bacterial swarming motility. 
Transcription of rsmZ is subject to the direct control of 
the GacS/GacA two-component system, and influenced by 
the sensor kinases RetS and LadS (21,22). Among them, it 
is the response regulator GacA that activates the rsmZ 
transcription by directly binding to the rsmZ promoter 
(17,23). Additionally, the H-NS-like protein MvaT acts 
as an rsmZ transcriptional repressor by binding to its 
promoter (23,24). Adding further to the complexity of 
the regulatory mechanisms that control motility and 
biofilm formation, evidence suggests the presence of 
some yet to be identified factors associated with the 
MvaT modulation of rsmZ expression (24). 

Given the central role of RsmA/RsmZ system in modu- 
lation of P. aeruginosa cell motility and biofilm formation, 
it is of considerable interest to investigate the detailed 
regulatory mechanisms that govern its expression. Here 
we identified a novel transcription factor BswR, which is 
involved in the regulation of P. aeruginosa swarming 
motility and biofilm formation. Overexpression of bswR 
enhanced the transcription of rsmZ, downregulated the 
expression of flagellar and T4P genes and attenuated the 
bacterial swarming motility. We also found that BswR 
binds to the promoter of rsmZ and counteracts the repres- 
sion of MvaT. Furthermore, we determined the crystal 
structure of BswR, which reveals that BswR forms a 
homodimer as the functional unit and provides clues on 
how BswR could interact with target promoter to regulate 
gene expression. These genetic and structural findings 
provide a new insight into the complicated and 
sophisticated regulatory mechanisms that govern the 
RsmZ/RsmA regulatory system and bacterial motility 
and biofilm formation. 



MATERIALS AND METHODS 

Bacterial strains, plasmids and culture conditions 

Bacterial strains and plasmids used in this study were 
listed in Supplementary Table SI. Both Escherichia coli 
and P. aeruginosa strains were maintained in Luria- 
Bertani (LB) broth with shaking at 250 rpm or on LB 
agar plates at 37°C. When necessary, antibiotics were 
included in medium as follows: 150 and 300mg-ml _1 car- 
benicillin, 5 and 50mg-ml _1 gentamicin and 5 and 
100mg-ml _1 tetracycline for E. coil and P. aeruginosa, 
respectively. 

Expression plasmid constructs were generated by 
standard methods and verified by DNA sequencing. The 
plasmids were transformed into E. coli by heat shock and 
P. aeruginosa strain by electroporation unless otherwise 



stated. To generate the fow7?-overexpressing construct 
pl9-bswR, the encoding region of bswR was amplified 
with polymerase chain reaction (PCR) primers 5 -GCGA 
ATTCTGCTAGGTACCCGGCTAAAAG-3 / and 5-GC 
GAATTCTCACAGTTCACTCCTTGTGC-3 / , digested 
with EcoRI and inserted into the corresponding site of 
pUCP19. The pl9-PA2781 was constructed similarly by 
using the PCR primers 5 -GCGAATTCTGAACAAAA 
GCTTGGATTCAG-3' and 5'-GCGAATTCCTCAGCC 
TCGCGCCAG-3'. The rsmZ reporter construct 
pME-rsmZ was generated by cloning the EcoRI-digested 
rsmZ promoter into the enzyme site of pME6016 (25). 
PCR primers for pME-rsmZ construction were 5'-GGA 
ATTCCAGTGACGGAAAACCTTAG-3' and 5 -GGAA 
TTCCCTGTACGCAGGAGTGATA-3'. To delete the 
16th-99th amino acids of BswR, the construct 
pEX18G-bswRl was prepared by overlapping PCR to 
generate the bswR allele bswRl, which lacks an internal 
fragment corresponding to nucleotides from 48th through 
297th base pair (bp). The mutated allele was 
integrated into the chromosome of P. aeruginosa by hom- 
ologous recombination, and the plasmid vector fragment 
was removed as previously described (26). The 
resultant mutants were verified by PCR and DNA 
sequencing. 

P-Galactosidase assays 

P-Galactosidase activity was measured in duplicates and 
repeated at least twice. Bacterial start cultures were diluted 
(1:1000) in LB, grown for 8h at 37°C and then aliquots 
were taken for measurement of enzyme activity. The 
P-galactosidase activity was quantified as previously 
described (27). 

Transposon mutagenesis 

The mariner transposon carried by plasmid pBT20 was 
used for mutagenesis of P. aeruginosa PAOl following 
the procedures as described (28). Mutants were screened 
on agar plates, which is a basic minimal nutrient medium 
supplemented with 0.2% mannitol as sole carbon source, 
0.2% ammonia sulphate as sole nitrogen source and gen- 
tamicin (50mg-ml _1 ) for selection of transposon mutants. 
The colonies showing smaller size than wild-type control 
were selected for further analysis. Arbitrary PCR, using 
primers of 5 / -GTCGASWGANAWGNA-3 / with 5'-GTG 
CAAGCAGATACGGT GACGAT-3' and 5 -TGACGA 
TCCCGCAGTGGCTCTC-3', was used to identify the 
genes disrupted by transposon insertion. The DNA se- 
quences flanking the transposon insertion were analyzed 
with the NCBI BLAST server (http://www.ncbi.nlm.nih. 
gov/blast/) and the Pseudomonas genome database (http:// 
www.pseudomonas.com/). 

Transcriptional profiling 

PAOl, PAOl (pi 9), AbswR and AbswR(bswR) cultures 
were grown overnight in LB medium. Bacterial cells 
were then diluted to OD 60 o of 0.05 and grown at 37°C 
with shaking at 250 rpm. When OD 600 reached ~1.0, 
cells were harvested and total RNAs were extracted with 
RNeasy Protect Bacteria Mini kit (Qiagen). DNA 
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contaminants were removed by DNase I. The quality of 
purified RNAs was analyzed by electrophoresis and 
quantified by Nanodrop 1000 Spectrophotometer 
(NanoDrop Technologies). Reverse transcription, frag- 
mentation and complementary DNA labeling were con- 
ducted as described by the manufacturer (Affymetrix). 
The processed samples were hybridized to the 
Affymetrix GeneChip of P. aeruginosa, and chips were 
washed, scanned and analyzed following the instructions 
from the manufacturer. Hybridization signals were pro- 
cessed using the statistics software MAS-5.0 from 
Affymetrix. 

Quantitative real-time reverse transcription polymerase 
chain reaction 

RNA samples were prepared as described above for 
transcriptional profiling analysis. An aliquot of 0.1 ug 
total RNAs was used as template, and real-time reverse 
transcription-PCR analysis was performed on Lightcycler 
vl.5 (Roche), using the SYBR® Green RT-PCR kit 
(Qiagen). The housekeeping gene rpoC, which encodes 
the (3 ; subunit of RNA polymerase, was included as an 
internal control. Results were presented as the ratio of 
target gene expression level versus the control gene (29). 
Real-time RT-PCR primers specific to rsmZ were 5'-CGT 
ACAGGGAACACGCAAC-3 / and 5 -ATTACCCCGCC 
CACTCTTC-3 / , and the primers for rpoC were 5-CTGTT 
CAAGCC GTTCATTTTC-3' and S'-CTTGATGGTGG 
TGGCCATA-3', respectively. 

Protein expression and preparation 

For biochemical assays, protein expression constructs were 
generated using the vector pET28a. For preparation of 
BswR, the encoding region of bswR was amplified with 
PCR primers y-CGGGATCCCTAGGTACCCGGCTAA 
AAGC-3' and 5 / -CGGATGCGTGTCGCGTAG-3 / , 
digested with BamHI and Hindlll and cloned into the cor- 
responding sites of pET28a (Novagen). For preparation of 
SpdH, the full-length spdH was amplified from the genomic 
DNA using the PCR primer pair 5'-CGGGATCCATGAC 
CATCTCCCGCCGCGA-3' and 5 -CCCAAGCTTCCGC 
TGCCGTCGCCGTCCT-3', and cloned into pET28a in 
the same way as described for bswR. After verification by 
DNA sequencing, pET28a-bswR and pET28a-spdH were, 
respectively, transferred into E. coli BL21(DE3), and the 
transformants were grown in LB media containing 
100mg-ml _1 kanamycin. Protein expression was induced 
by addition of 0.1 mM IPTG at 18°C for overnight. Cells 
were then harvested by centrifugation (4500 rpm, 30min, 
4°C), disrupted by French pressure cell, and the lysates 
were cleared by ultracentrifugation (25 000 rpm, 1 h, 4°C). 
For preparation of MvaT, the full-length mvaT was first 
cloned into pET28a using PCR primers 5'-GCGCATATG 
TCCCTGATCAACGAATATCG-3' and 5'-GGGCTCGA 
GGCCGAGCAGGGTGGCCCAG-3', and then the 
fragment carrying mvaT and His-tag was inserted into the 
pUCP19 vector to obtain pl9-mvaT. After sequence verifi- 
cation, the vector pl9-mvaT was then transferred into 
AmvaT(P r5mZ -lacZ) for protein purification. His 6 -tagged 
MvaT was then purified from the P. aeruginosa strain 



using Ni -affinity chromatography as recommended by 
Qiagen, followed by dialysis against the electrophoretic 
mobility shift assay (EMSA)-binding buffer (10 mM Tris, 
pH 7.5, 50 mM KC1, 1 mM dithiothreitol). Protein purity 
and integrity were examined by sodium dodecyl sulphate- 
polyacrylamide gel electrophoresis, and the concentration 
was determined by Bradford protein assay (Bio-Rad). 

For crystallographic analysis, bswR was amplified by 
PCR and cloned into modified pET26b(+) vector 
(Novagen), in which BswR was expressed as a 
C-terminal 6xHis-tagged fusion protein. The expression 
vector pET26M-bswR was then transformed into E. coli 
BL21(DE3) and expressed as described above. The BswR 
protein was first isolated by Ni-NTA column (GE 
healthcare) under a linear gradient elution of 0.0-0.5 M 
imidazole in the buffer 50 mM NaH 2 P0 4 , pH 7.5, 
300 mM NaCl, 5mM (3-mercaptoethanol ((3-ME), and 
further purified by HiLoad superdex 75 26/60 gel filtration 
column (GE healthcare), with the elution buffer 20 mM 
Tris, pH 7.5, 200 mM NaCl and 5mM [3-ME. The 
purified proteins were concentrated to 5.2mg«ml _1 and 
stocked in the buffer (20 mM Tris pH 7.5, 200 mM 
NaCl, 5mM (3-ME) for crystallization trials. 

Electrophoretic mobility shift assay 

The promoter region of rsmZ was labeled with biotin by 
PCR using primers S'-Biotin-TTCCAGTGACGGAAAA 
CCTTAG-3' and S'-Biotin-TTCCCTGTACGCAGGAG 
TGATA-3'. The control probe pexsC that extended the 
promoter region of exsC from P. aeruginosa PAOl was 
PCR-amplified using S'-biotinylated primers 5'- CTCCTA 
AAGCTCAGCGCATG-3' and 5 ; -GTCAGTCCTATTT 
CACCCAGAG-3'. The biotinylated promoter was 
purified with QIAquick PCR Purification kit (Qiagen) 
and quantified by Nanodrop spectrophotometer. EMSA 
was carried out using LightShift® Chemiluminescent 
EMSA kit following the manufacturer's instruction 
with modifications (Pierce). Each reaction (10 ul) 
contained the biotin-labeled probe (0.2 ng), 1 jig poly 
(dI»dC) (2'-deoxyinosinic-2'-deoxycytidylic acid) and 
various amounts of His 6 -BswR protein as indicated. In 
parallel, His 6 -SpdH was included in the assay as control. 
The reaction mixtures were incubated at 25°C for 20min, 
and then the samples were separated on 6% polyacryl- 
amide gels. The DNA probes were then transferred onto 
a nylon membrane, cross-linked by ultraviolet light at 
254 nm and detected by chemiluminescence as described 
by the manufacturer (Pierce). 

Motility assay 

The motility of P. aeruginosa was examined on LB agar 
(1.5%) plates at 37°C. Unless stated otherwise, 2|il of 
fresh P. aeruginosa cultures was spotted in the center of 
plates and grown at 37° C for 2 days. Swarming, swimming 
and twitching motilities were respectively determined as 
described previously (4). Results were recorded after 
incubation overnight at 30° C. 
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Biofilm assay 

Biofilm formation was determined in 15-ml borosilicate 
tubes. Briefly, 1 ml of LB medium was inoculated to a 
final OD 6 oo of 0.1 and statically incubated at 37°C for 
6h. Biofilm cell mass was virtualized by staining with 
0.1% crystal violet and quantified by measuring the ab- 
sorbance at 595 nm after solubilizing crystal violet in 96% 
ethanol. Each experiment was independently repeated for 
at least three times. 



DNase I protection assay 

The promoter region of rsmZ was amplified by PCR with 
FAM-labeled primers 5 / -(6-FAM)-ATTCTGGAGAAGA 
ATGGCCTGTC-3 / and S'-TTGCGTGTTCCCTGTACG 
C-3' from the genomic DNA of P. aerugionsa PAOl. The 
experiment was done as described previously (30) in four 
repetitions with some modifications. Briefly, 100 ng of 
DNA was diluted in 20 ul of gel shift buffer, then 5 jil of 
purified BswR and 1 3 (il of water were added to the 
reaction; binding was allowed to proceed for 20min at 
room temperature. Following several trials, 10U of 
DNase I (Roche) was added and incubated for 5min, 
and then the enzyme was heat-inactivated for 5min at 
95°C; DNA was purified by using Qiagen PCR kit and 
eluted on 25 ul of water for further analysis. 
Dideoxynucleotide-based sequencing was carried out 
using Thermo Sequenase Dye Primer Manual Cycle 
Sequencing kit (USB) according to the manufacturer's in- 
structions, and samples were analyzed by using the 
Applied BioSystems 3730x1 DNA analyzer. 



Crystallization and structure determination 

Crystallization trials were carried out using sitting-drop 
vapor diffusion method by Robot Phoenix (Art Robbins 
Instruments) at 20° C. The diffraction quality crystals of 
BswR were obtained in the drop by mixing protein sample 
and reservoir solution (0.1 M HEPES, pH 7.5, 20% w/v 
polyethylene glycol 10 000) with a volume ratio of lOOnl 
to 200 nl. For data collection, crystals were cryoprotected 
in reservoir solution supplemented with 20% glycerol, and 
flash-frozen in liquid nitrogen. 

Diffraction data set was collected in Bruker AXS X8 
Proteum X-ray system (wavelength = 1.54178 A) (Bruker 
Axs Inc., USA). The data set was indexed, scaled and 
integrated with iMosflm (31). Using 3LFP [Protein Data 
Bank (PDB) code] as a search model, initial phases were 
obtained by molecular replacement using program Balbes 
in CCP4 suite (32), followed by model building with ARP/ 
wARP (33). Model visualization and manual model 
building were carried out with COOT (34), and the refine- 
ment was done with CNS (35). The final model was 
examined using PROCHECK (36), crystallographic data 
and refinement statistics are listed in Table 1. Protein 
interface areas were calculated using PISA sever (37). 
Structural illustrations were generated using PyMol 
(DeLano Scientific). 



RESULTS 

BswR regulates P. aeruginosa motility and biofilm 
formation 

To identify unknown regulatory factors associated with 
the motility of P. aeruginosa, a transposon mutant 
library was screened for alteration in colony size on LB 
agar plates. Of ~2000 mutants screened, three mutants 
with colony size smaller than the wild-type control were 
selected and designated as Motl, Mot2 and Mot3 
(Supplementary Figure SI A). Sequence analysis of the 
three mutants showed that transposon was inserted in 
PA2285 (Motl), PA3966 (Mot2) and the intergenic 
region of PA2779 and PA2780 (Mot3), respectively. 
Among them, PA2285 encodes the transcriptional factor 
AmrZ associated with the regulation of P. aeruginosa 
motility (38), and PA3966 encodes a hypothetic protein 
with unknown function, whereas PA2779 and PA2780 
encode a hypothetic protein and a putative 
DNA-binding protein, respectively. The mutant Mot3, in 
which the transcriptional expression of PA2780 was 
significantly upregulated by transposon insertion 
(Supplementary Figure S2), was selected for further 
analysis in this study. 

As shown in Figure 1A, the ORFs PA2779 and PA2780 
were divergently orientated, sharing an intergenic region 
of 407 bp as putative promoters. The mariner transposon 
was inserted at the 12th bp upstream of PA2780, and the 
transposon-borne tac promoter was orientated in the same 
direction as that of PA2780, which is known to increase 
the transcriptional expression of downstream genes (28). 
Following PA2780 is PA2781, a 339-bp ORF, which 
overlaps the last 1 bp of PA2780 and possibly forms an 
operon with PA2780. Similar to PA2780, PA2781 also 
encodes a hypothetic protein with unknown function. 
No discernible motifs were found in both proteins, but 
bioinformatics analysis suggested that PA2780, 
comprising 114 amino acid residues, possesses a potential 
xenobiotic response element (XRE)-type DNA-binding 
domain at its N- terminus (Figure 1A). According to its 
role in bacterial swarming regulation as described below, 
PA2780 was designated hereafter as bswR (Bacterial 
SWarming Regulator). 

In other bacterial species, the proteins containing the 
XRE-type DNA-binding domain have been implicated 
in regulating plasmid copy number and transcriptional 
expression of bacteriophage genes (39). However, the 
biological function of BswR in P. aeruginosa has not yet 
been characterized. To elucidate its role in bacterial 
motility, we generated the mutant AbswR by in-frame 
deletion of bswR in wild-type strain PAOl. We first 
checked the colony size on LB agar plates and found the 
deletion mutant formed much larger colonies than the 
wild-type and the fow7?-overexpressing derivative 
(Supplementary Figure SIB), which was not due to 
altered bacterial growth rate (Supplementary Figure S3). 
We then examined the bacterial motility on specific agar 
plates for analysis of swarming, swimming and twitching 
motility (4). As shown in Figure IB, we found deletion of 
bswR dramatically increased the bacterial swarming 
motility and the changed phenotype was rescued by 
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Table 1. Summary of crystallographic data and refinement statistics 
Data collection 



Space group 


C2 


Unit cell dimensions 




a,b,c (A) 


a = 63.0, b = 43.9, c = 35.0 


a,P,y (°) 


a = 90, P = 115.2, y = 90 


Resolution (A) 


30.0-2.3 (2.42-2.30) 


Rsym (%) 


8.5 (22.1) 


I/ctI 


6.4 (2.9) 


Completeness (%) 


97.3 (83.5) 


Refinement 




T? f^cr^lnti r\r\ ( W 
rVCoVJlLALUJll y-f^/ 


30.0-2.3 


ATnmnpr of nnimip rpflpptinnc 


3818 


K work/ K free v /o ) 


11 Q/98 f\ 






Protpm 


803 


Water 


50 


Average B factor (A 2 ) 




Main chain 


26.8 


Side chain 


29.9 


Waters 


41.3 


Rmsd 




Bond length (A) 


0.005 


Bond angle (°) 


0.978 


Ramachandran plot 




Most favored regions (%) 


94.5 


Additional allowed regions (%) 


4.4 


Generously allowed regions (%) 


1.1 


Disallowed regions (%) 


0.0 



and <I> is the mean intensity for that reflection. 

Reflections with I>a was used in the refinement. 

c R WO rk = I Fobs - F ca i c |/|F obs | where F calc and F obs are the calculated and 

observed structure factor amplitudes, respectively. 

d R free = as for R wor k, but for 9.4% of the total reflections chosen at 

random and omitted from refinement. 

individual B-factor refinements were calculated. 



in trans expression of the wild- type bswR but not PA2781 
in the deletion mutant AbswR (Figure IB). We also tested 
the effect of bswR overexpression on bacterial swimming 
and twitching motility, but no significant change was 
observed (data not shown). These data suggest that 
BswR is specifically associated with the regulation of 
P. aeruginosa swarming motility. 

Given that bacterial motility is critical for biofilm 
formation in P. aeruginosa (40), we tested the impact of 
bswR on biofilm formation. As depicted in Figure 1C, 
biofilm formation was almost completely abolished in 
the deletion mutant AbswR. In contrast, overexpression 
of BswR in the deletion mutant AbswR dramatically 
increased the biofilm formation. Surprisingly, overexpres- 
sion of PA2781 in AbswR had no effect on bacterial 
swarming motility (Figure IB), whereas this resultant 
strain produced more biofilm than that produced by the 
wild-type PAOl (Figure 1C), suggesting that the hypo- 
thetical protein PA2781 may play a specific role in 
biofilm development, which is currently under 
investigation. 

BswR modulates the biosynthesis of flagella and 
type IV pili 

To understand how BswR could regulate the swarming 
motility, we conducted microarray analysis to compare 



< pa2779 T 



I bswR (pa278dT S X pa2781 > I po2782 > 



HTH 

XRE 



PAOl 



Mot3 



AbswR(bswR) 



AbswR 



AbswR(pa2781) 



3.0 



Q 2.5 



2.0 



■2 1.0 



0.5 



0.0 



UJi 



Figure 1. BswR modulates the motility and biofilm formation of 
P. aeruginosa. (A) Genetic organization of bswR. Arrow of ORF 
indicates the transcriptional direction. The solid diamond indicates 
the location of transposon insertion in the mutant Mot3 and the 
transcription direction of the transposon-borne F tac promoter is 
indicated by a solid arrow. The domain structure of BswR was 
predicted by using the SMART program (http://smart.embl-heidel 
berg.de/). (B) Representative swarming motility image of 
P. aeruginosa wild-type strain PAOl and its derivatives. (C) Biofilm 
formation of PAOl and its derivatives. 



the global gene expression profiles of the deletion 
mutant AbswR and wild-type strain PAOl. The results 
showed that deletion of bswR significantly changed the 
expression patterns of 454 genes (>2-fold), with 258 
being upregulated and 196 downregulated (Supplementary 
Table S2). Considering the critical roles of flagella and 
T4P in bacterial motility and biofilm formation, we spe- 
cifically compared the expression patterns of the genes 
associated with flagella and T4P in PAOl and AbswR, 
and found that >40 flagella and T4P genes, including 
the gene clusters flgB-flil, cheV-flgN, fimT-pilE and pilQ- 
pilN, were upregulated (>2-fold) in the deletion mutant 
(Figure 2 A and B). These upregulation patterns were 
verified by real-time RT-PCR analysis (Supplementary 
Figure S4). Consistently, overexpression of bswR in the 
strain AbswR resulted in decreased (>2-fold) transcrip- 
tional expression of 14 flagella and 15 T4P genes 
compared with the wild-type PAOl (Supplementary 
Figures S2 and S5). These results suggest that BswR regu- 
lates swarming motility likely through modulation of 
flagella and T4P biogenesis. 
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pilH pilA pilB pilD pilR fimT fimU pilV pilW pilYl pilQ pilO pilN 



C pelA pelB peIC pelD pelF peIG 




Figure 2. BswR regulates the transcription of the genes associated with the biogenesis of flagella (A), type IV pili (T4P) (B) and biofilm (C). 
Microarray analyses were performed using the exponential phase bacterial cells grown in LB medium. The BswR-regulated genes (>2-fold) 
involved in the formation of flagella, T4P and biofilm were presented and organized according to their corresponding operons with arrows indicating 
the transcriptional direction. The transcriptional fold changes were presented for the genes increased in the mutant AbswR in comparison with 
wild-type strain PAOl. 



In addition to flagella and T4P, extracellular polysac- 
charide (EPS) is another factor associated with biofilm 
formation in P. aeruginosa (41); we then analyzed the 
genes responsible for EPS production based on micro- 
array data. Among the two gene clusters, i.e. the pel 
operon and the psl operon, known to be associated with 
EPS production, we found that only the expression of pel 
operon was significantly decreased in the AbswR mutant 
compared with the wild- type strain (Figure 2C). These 
data are in line with the previous study that mutation 
of the pel genes results in decreased biofilm formation 
(42) and establish a link between BswR and the pel 
genes in biofilm formation. Taken together, these results 
suggest that the influence of BswR on bacterial biofilm 
formation is not only related to its negative role in regu- 
lation of flagella and T4P biogenesis, but may also be 
associated with its positive effect on modulation of EPS 
production. 



BswR regulates the transcription of small RNA rsmZ 

Besides the genes related to flagella, T4P and EPS 
biosynthesis, genes involved in protein secretion were 
also regulated by bswR (Supplementary Table S2). As 
shown in Figure 3A, expression of the genes for type III 
secretion system (T3SS) (including exoT, exsE, pscRQPO, 
per 1234, popB and pscEF) was dramatically increased by 
deletion of bswR. Intriguingly, we noted that these 
expression patterns were mutually reciprocal with those 
caused by the deletion of the RNA-binding protein 
encoded by gene rsmA (20), and reminiscent to those 
caused by double deletion of small RNAs rsmZ/Y 
(Figure 3 A) (23). It has been reported previously that 
deletion of rsmZ/Y results in increasing expression of 
T3SS genes and decreasing expression of T6SS genes, 
whereas deletion of rsmA increases the expression of 
T6SS genes but decreases the expression of T3SS and 
T4P genes (20,23). 
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Figure 3. BswR regulates the transcription of small RNA rsmZ. (A) The gene expression patterns caused by bswR deletion were similar to those 
caused by rsmZ/Y deletion. Y-axis: the fold-changes of T3SS genes in AbswR compared with the wild-type and those in ArsmZ/Y as published 
previously (23). X-axis: Representatives of the TTSS genes whose mRNA levels were altered by deletion of bswR and rsmY/Z in P. aeruginosa, 
respectively. (B) Deletion of bswR resulted in decreased expression of rsmZ. The P-galactosidase activity conferred by FrsmZ-lacZ promoter fusion 
was determined when bacterial cells grown in LB medium reached at stationary phase. (C) The relative mRNA level of rsmZ revealed by real-time 
qRT-PCR analysis. 



Given that the RsmZ/7-RsmA system is known to be 
associated with the regulation of P. aeruginosa motility 
(18), the above results motivated us to postulate that 
BswR might regulate the expression of rsmZ/Y. As no 
specific probes matching rsmZ/Y were designed in the 
GeneChips from Affimatrix, we tested this hypothesis by 
generating two reporter constructs by fusing the pro- 
moters of rsmZ (Prsmz) and rsmY (F rsmY ) separately to 
the coding region of lacZ and introduced these constructs 
into PAOl and its derivatives. Analysis of the 
Prsmz-directed P-galactosidase activity showed that 
deletion of bswR decreased the rsmZ expression, and 
overexpression of bswR increased the rsmZ expression 
level (Figure 3B). However, such differences were not 
observed when the P r5my -lacZ fusion gene was assayed 
(data not shown). These results suggest that BswR 
regulates the transcriptional expression of rsmZ but not 
rsmY. Consistent with these results, real-time RT-PCR 
assay also showed that the transcript level of rsmZ was 
decreased by ~60% in the fow7?-deleted mutant and 
increased by ~ 3-fold in the Z?swi?-overexpression strain 
(Figure 3C). 

BswR requires GacA and MvaT for upregulation of rsmZ 

GacA, a response regulator of the two-component system 
GasA/GacS in P. aeruginosa, has been well known to 



regulate rsmZ by directly binding to its promoter. It has 
been shown that disruption of GacA by transposon 
significantly reduced the rsmZ expression (23). To charac- 
terize the BswR-mediated regulation of rsmZ, we 
compared the transcription level of rsmZ in the wild-type 
strain and the relevant mutants. Consistent with previous 
report (23), disruption of gacA substantially reduced the 
transcriptional expression of rsmZ compared with the 
wild-type strain (Figure 4A). Notably, while in trans 
expression of bswR in wild-type strain increased the 
transcriptional expression of rsmZ, expression of bswR 
in the gacA mutant strain had no effect on rsmZ 
transcription (Figure 4A). The findings were further 
consolidated by examination of phenotype changes. As 
depicted in Figure 4B, overexpression of bswR inhibited 
the swarming motility of the wild-type PAOl. In contrast, 
the influence of overexpressed BswR was diminished when 
gacA was mutated, suggesting that BswR requires a 
functional GacA for the upregulation of rsmZ 
transcription. 

In addition to GacA, the H-NS-like DNA-binding 
protein MvaT also binds to the rsmZ promoter and 
negatively regulates the transcription of rsmZ (24). To 
test whether BswR requires MvaT for the upregulation 
of rsmZ, we compared the transcription level of rsmZ in 
the wild- type and the relevant mvaT mutant strains. 
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Figure 5. BswR binds to the promoter region of RsmZ. (A) EMSA of 
BswR binding to the promoter of rsmZ. EMSA assays were performed 
in presence of various amounts of His6-tagged BswR using the 
biotinylated promoter of rsmZ (P-rsmZ) as described in the text. 
His 6 -tagged SpdH, a protein associated with spermidine metabolism, 
and the biotinylated promoter of exsC (PexsC) from P. aeruginosa 
PAOl, were included as controls. (B) Identification of BswR binding 
sites in the promoter region of rsmZ using Dye Primer Sequencing on 
an automated capillary DNA analyzer. Top: electropherograms 
showing the protection region (red line) of the rsmZ promoter after 
digestion with DNase I following incubation in the presence (BswR) or 
absence (CK) of BswR. (C) Sequence analysis of the promoter region 
of rsmZ. Sequences for PCR primers were italicized. Binding sites for 
BswR as determined in this study were underlined in red; binding sites 
for GacA and MvaT were underlined in blue and green, respectively. 
The arrow indicated the transcription start site of rsmZ. 



Consistent with the previous study (23), deletion of mvaT 
resulted in > 2-fold increase in rsmZ transcription 
compared with the wild- type strain PAOl (Figure 4A). 
In contrast, overexpression or deletion of the bswR gene 
in the background of mv^T-deletion mutant did not cause 
any change at the transcript level of rsmZ (Figure 4A). 
As a control, overexpression of bswR in PAOl inhibited 
the swarming motility of the wild-type strain, and this 
inhibition was almost completely relieved in the mvaT 
mutant (Figure 4B). Given that overexpression of bswR 
increases the rsmZ transcription in the wild-type strain, 
these results suggest that BswR may act by counteracting 
the repressor MvaT in upregulation of rsmZ. Consistent 
with this notion, we found that the double deletion mutant 
AbswRAmvaT had a comparable transcription level of 
rsmZ as that of the single deletion mutant AmvaT 
(Figure 4A). In line with these findings of genetic 
analysis, the swarming assay showed deletion of bswR in 
strain PAOl increased the bacterial motility, but such an 
effect was vanished when mvaT was also deleted (Figure 
4B). 



BswR binds to the promoter region of rsmZ 

To characterize the function of BswR, we examined 
whether BswR could directly interact with the promoter 
of rsmZ in vitro. First, we expressed and purified the 
recombinant BswR as a His 6 -tagged protein and 
conducted the EMSA. As shown in Figure 5A, the 
mobility shift of the rsmZ promoter was affected by 
BswR in a dosage-dependent manner, demonstrating 
BswR binds to the rsmZ promoter. In contrast, the 
control protein SpdH, which is a spermindine 
hydrogenase in P. aerugionsa (43), did not interact with 
the rsmZ promoter. Furthermore, the purified BswR did 
not bind to the control DNA fragment PexsC, which is the 
promoter of gene exsC in P. aeruginosa (Figure 5A). 
Cumulatively, these results indicate that BswR positively 
regulates the transcriptional expression of the small RNA 
rsmZ by specifically binding to the promoter region of 
rsmZ. 

To identify the binding sites of BswR in the rsmZ 
promoter region, we conducted the DNase I protection 
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assay by using the PCR-amplified promoter region of 
rsmZ, which is labeled with the fluorescent dye FAM 
(fluorescein amidite). The results revealed that the DNA 
stretch 'ACCCGGCGTTCGCCGGGT' was clearly 
protected from the DNase I digestion on incubation 
together with BswR (Figure 5B). This protected region 
appeared as a palindrome, typical of the transcriptional 
factor binding sites in bacteria. In the promoter region of 
rsmZ, it extended from positions —226 to —208 relative to 
the transcriptional start site of rsmZ, upstream of the 
binding regions of GacA and MvaT (Figure 5C). 
Moreover, sequence analysis showed this protected 
region is absent in the promoter region of rsmY. This 
result was consistent with our findings that BswR has 
little effect on the transcription of rsmY, further 
supporting that BswR specifically regulates rsmZ by 
directly binding to its promoter region. 

Overall structure of monomeric BswR 

To further understand the molecular mechanisms with 
which BswR interacts with the promoter of rsmZ, we 
proceeded with structural analysis by X-ray crystallog- 
raphy. The crystal o structure of BswR was determined 
and refined to 2.3 A resolution. While 101 residues were 
visible in the final model, the last 13 residues and the 
6xHis-tag in C-terminus were not well ordered in the 
electron density map and thus excluded. Then, the BswR 
protein model was refined and water molecules were 
located, finally giving R-work/R-free factor to 22.9/ 
28.6%, respectively. The final model was examined using 
PROCHECK: 94.5% of the residues lie in the most 
favored regions, with the remaining 4.4 and 1.1% 
assigned to the additionally and generously allowed 
regions, respectively. The final refinement statistics for 
the structure are listed in Table 1 . 

The overall topology of BswR monomer contains 
seven helices, of which five form a compact helix bundle 
and the other two assemble as extended C-terminal helices 
(Figure 6A). Consistent with the prediction, our structure 
reveals that the N-terminus of BswR comprises 
oc2-turn-oc3, a typical helix-turn-helix (HTH) motif for 
DNA-binding (in which oc2 is for stabilization and oc3 
for recognition). The two helices (oc6 and oc7) with a 
V-shaped arrangement in the C-terminus vary largely 
across different organisms, which meets the requirement 
for diverse ligand binding. 

Dimerization of BswR 

Based on the gel filtration profile during protein 
purification, we found that BswR exists as a dimer in 
solution. This was further confirmed from the crystal 
structure: the two monomers with the V-shaped 
arrangement of the C-terminal helices intertwine together 
to form a homodimer by crystallographic symmetry. Using 
PISA and Contact (CCP4 suite) programs (33,39), we 
investigated the dimerization interface of BswR. In either 
monomer, 94 of 101 amino acids were exposed on the 
surface. On dimerization, an average interface area of 
1438.7 A 2 was buried, accounting for 20.8% of the 
solvent-exposed surface area (6923.4 A 2 ) of each 



monomer. In total, 45 contacts were identified between 
two monomers with an intermolecular distance of 
<3.5A, the majority of which were located at the 
C-terminus. Furthermore, dimerization was mainly 
mediated through hydrophobic interaction, which 
involves the residues Leu71, Leu74, Leu81, Phe85, Leu89, 
Thr90, Ile93, Leu97 and Ala99 (Figure 6D and E). In 
addition, the intermolecular hydrogen bonds between the 
main chain N atom of Ala62 and the main chain O atom of 
Tyr65 (Ala62-Tyr65 / , Tyr65-Ala62 / , where the primes refer 
to the second monomer, same as below), the salt bridges 
(Lys77-Glu94, Glu94-Lys77 / ), as well as the aromatic k-k 
stacking interactions (Tyr65-Tyr65 ; and Phe92-Phe92 / ) 
(Figure 6C and D), were found to contribute to the BswR 
dimerization. In the dimeric structure, BswR exhibits an 
uncommon large area of dimer interface, which is formed 
by the extension of C-terminal helical region (helices 5-7, 
residues 61-101). Clearly, the dimerization of BswR was 
mainly achieved by the interleaving of the two extended 
helices in its C-terminus, although the parallel helix 5 also 
contributes to the dimer formation. 

XRE comprises the second-most frequently occurring 
regulator family in bacteria that control diverse biological 
functions, and these regulators share a characteristic 
N-terminal DNA-binding domain with a HTH motif, 
whereas the C-terminal region is extremely variable. 
Structural comparison of BswR with other XRE family 
proteins, restriction-modification (R-M) system control- 
ler proteins (C protein) Csp231I from Citrobacter sp. 
RFL231 (44), SinR from Bacillus subtilis (45), C protein 
Espl396I from Enterobacter sp. RFL1396 (46) and EspR 
from Mycobacterium tuberculosis (47), shows that they 
share a similar fold for DNA-binding domain 
(Supplementary Figure S6), despite the variations mainly 
in the turn of HTH (oc2-turn-oc3) which could be of 
importance for orientating the recognition helix oc3 for 
an optimal affinity for DNA (47). Except for C protein 
Esp 13961 with one C-terminal helix, these proteins contain 
two C-terminal extended helices forming a V-shaped ar- 
rangement, which is involved in dimer assembling. The 
orientation of this V-shaped arrangement relative to the 
similar N-domain is varied, which consequently impacts 
on its DNA-binding (47). In contrast to the structural 
similarity, the sequence identity is relatively low 
(Supplementary Figure S6). 

Structural basis of BswR for DNA binding 

BswR possesses a calculated isoelectric point (pi) of 9.2, 
indicating that it is positively charged under physiological 
condition. EMSA assay showed that BswR can bind to the 
promoter region of rsmZ, and DNase I protection assay 
further revealed that ' ACCCGGCG TT CGCCGGGT ', a 
palindromic sequence (underlined) in the rsmZ promoter, 
is the binding site of BswR. To elucidate how BswR binds 
to DNA, a structural alignment of BswR onto Esp 13961- 
DNA (PDB ID: 3CLC) resulted in a DNA-binding model 
(Figure 7A), in which no steric clash between BswR and 
the double-stranded DNA (dsDNA) was observed. In this 
model, the HTH motif in each monomer of BswR dimer is 
the major component to interact with the dsDNA, and the 
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Figure 6. The structure of BswR and its dimerization. (A) Topology of the BswR monomer. The N- to C-termini of BswR is colored from blue to 
red. (B) The structure of BswR dimer. In a homodimer, one monomer is colored green, and the second brown (where the primes referred to the 
second monomer, same as below). The dimerization interface is framed by the black box, depicted as line, dashed line, and dotted line, respectively. 

(C) Enlarged view of dimer interaction along the parallel oc5. The residues involved in dimerization interface are shown as stick model with labels. 
Hydrogen bonds between A62-Y65' and A62'-Y65 were indicated as dashed lines, %- % stacking Y65-Y65' was indicated with a double-ended arrow. 

(D) Enlarged view of dimer interface along the V-shaped terminal. Residues involved in dimerization were labeled and shown as stick model. Salt 
bridges (K77-E94' and K77'-E94) are indicated as dashed lines. (E) Detailed interaction between a7 and al' . Residues involved in dimerization are 
shown as stick model, n- n stacking F92-F92' was indicated with a double-ended arrow. 



recognition helix oc3 was inserted into the consecutive 
major groove of the dsDNA, in a canonical way to estab- 
lish contact. 

We also mapped the electrostatic surface of BswR 
dimer, as shown in Figure 7B. The strong positive 
charge was mainly located at the N-terminal DNA 
binding domain, and the remaining surface charges were 
evenly distributed. The dsDNA molecules were laid across 
the positive-neutral-positive patches in BswR dimer, and 
most positively charged region colored blue was covered, 
thereby demonstrating that the model is nicely consistent 
with the surface charge distribution. Compared with other 
XRE family transcriptional factors, several conserved 
residues were mapped out for BswR to bind DNA 



(Figure 7C). These residues, including Arg34, Gln37, 
Arg42, His43, Asn46 and Arg48 in the N-terminal 
domain, face toward the dsDNA and probably contribute 
to the direct contact with DNA. 



DISCUSSION 

Bacterial motility and biofilm formation has been exten- 
sively investigated (3,8-11). In P. aeruginosa, swarming 
and biofilm formation are mutually exclusive, could be 
regulated by nitrogen limitation, controlled by quorum 
sensing and governed by a large number of transcriptional 
regulators (8,9). In this study, we discover a new layer of 
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Figure 7. A proposed model of BswR interacting with DNA. (A) Model depicting BswR-DNA complex. Structural alignment of one protein dimer 
of Espl396I-DNA complex (PDB ID: 3CLC) onto BswR dimer by fitting the two N-terminal domains was carried out with coot, subsequently the 
excess nucleotides far apart from BswR were removed. The dsDNA was shown as cartoon with sticks, same as (B). (B) A representative portion of 
the BswR-DNA complex structure with exposure of the HTH regions. Electrostatic potentials were rendered on the protein surface with a color scale 
between -40.0 kT e" 1 (Red) and +40.0 kT e" 1 (blue). The DNA was shown as cartoon with sticks. (C) The putative residues of BswR involved in 
DNA binding. The residues are predicted to interact with DNA are labeled and shown as sticks. For clarity, the DNA was shown as cartoon without 
sticks. 



regulatory cascade that tailors swarming and biofilm for- 
mation in P. aeruginosa, by identification of a previously 
uncharacterized protein BswR (Figure 1). Genome-wide 
transcriptomic analyses show that BswR specifically 
controls the activity of flagella and T4P (Figure 2). 
Alteration of the bswR expression affected the swarming 
of P. aeruginosa and also affected the expression of genes 
involving the assembly of flagella and the biogenesis of 
T4P, indicating that BswR regulates bacterial swarming 
likely through compromising the function of flagella and 
T4P. However, mutation of bswR did not significantly 
influence the flagella-dependent swimming and the 
T4P-dependent twitching, suggesting that BswR have 



additional roles in regulation of bacterial swarming 
under these experimental conditions. 

Our results provide new insights into how the small 
RNA rsmZ is regulated in response to cellular physiology 
and environmental stimuli (Supplementary Figure S7). 
RsmZ is a multifunctional regulator in P. aeruginosa. By 
counteracting the activity of RsmA, RsmZ acts as a 
negative regulator for bacterial swarming and a positive 
regulator for biofilm formation (18-20). In present study, 
we found overexpression of bswR inhibits bacterial 
swarming but favors biofilm formation, compatible with 
the previous results that overexpression of rsmZ resulted 
in loss of the swarming ability and increase of biofilm 



4574 Nucleic Acids Research, 2014, Vol. 42, No. 7 



production (23). Identification of rsmZ as the target of 
BswR provides a molecular mechanism through which 
BswR regulates the bacterial swarming and biofilm 
formation. 

Transcription of rsmZ is activated by GacA but 
repressed by MvaT, both physically binding to the 
promoter region of rsmZ (17,18,22,23,48). While the 
activity of GacA is titrated by a few two-component 
sensors, the repression of MvaT could be attenuated by 
some yet unidentified transcription factors (23). Our 
results provide a few lines of evidence suggesting that 
BswR is one of such factors. First, BswR activates the 
transcription of rsmZ (Figure 3). Second, BswR requires 
GacA for its activation (Figure 4). Third, this activation 
depends on the presence of MvaT (Figure 4). Fourth, 
BswR directly binds to the rsmZ promoter (Figure 5). 
Even so, the detailed mechanism through which BswR 
relieves the repression of MvaT requires further investi- 
gation. In the rsmZ promoter, binding site of GacA (TGT 
AAGCATTAACTTACA) is upstream of that of MvaT 
(AATTTAAGTAACTTATTGAAATTAAAGGAATAT 
ATTTTTCTAAAAT), separated by a DNA stretch with a 
length of 56 bp (Figure 5C). DNA footprint unravels that 
BswR binds to a separate site from that of MvaT, 
excluding the possibility that BswR and MvaT directly 
compete the binding sites for expression of RsmZ. In the 
crystal structure, BswR contains a positively charged 
DNA-recognition extension that dramatically deforms 
the promoter region during DNA binding (Figure 7), 
indicating that BswR may activate rsmZ by influencing 
the binding affinity of MvaT to the promoter region. 
However, previous results demonstrate that MvaT could 
associate with two small regulatory peptides for its repres- 
sion activity (49), and moreover, MvaT functioning as a 
member of H-NS family protein usually causes DNA 
bending and thus facilitates distant protein-protein inter- 
actions (50). It is possible that BswR counteracts MvaT by 
directly binding to each other. In addition, the MvaT of 
P. aeruginosa has been reported to form both dimers and 
higher-order oligomers and the binding of MvaT to DNA 
targets depends on the oligomerization of MvaT dimers 
(50). Therefore, it is also possible that BswR counteracts 
MvaT by interfering with the protein-protein interactions 
between MvaT dimers. 

Bacterial motility and biofilm are most important viru- 
lence determinants and survival strategies for microbial 
pathogens to cope with harsh environments and infecting 
hosts. Expression of these traits is under stringent 
regulation, and it responds to largely unidentified 
environmental signals. Although the regulatory cascade 
has been extensively studied, the environmental signals 
remain largely unknown. Structural comparison of 
BswR with other XRE family proteins reveals a similar 
fold for the N-terminal domain (Supplementary Figure 
S6). Further structural analyses suggest BswR functions 
as a homodimer, which uses positively charged helix to fit 
the major groove of the dsDNA for binding (Figure 7). In 
particular, BswR shares the most similar structure to C 
protein Csp231I, o with a root-mean-square deviation 
(rmsd) of 1.169 A for the main chain (Supplementary 
Figure S6). In Citrobacter sp. RFL231, Csp231I plays a 



key role in the temporal regulation of gene expression in 
bacterial R-M systems and are important mediators of 
horizontal gene transfer (44). The biological significance 
of this structural similarity is not known yet. Different 
from the N-terminal DNA-binding domain, the 
C-terminal region is structurally diversified with no 
typical homologs identified in the database. In the 
crystal structure, it comprises two alpha-helices with a 
V-shaped arrangement for dimerization. It is unknown 
whether this arrangement also involves the potential 
ligand binding for BswR, given that no ligand has been 
identified for these proteins, and little information is 
available on how the ligand-binding site could be 
formed. Analyses of previously published genome-wide 
microarray results show that the expression of BswR is 
not significantly influenced under the conditions of 
nutrient depletion, altered quorum sensing and increased 
temperatures, all of which have been implicated in 
bacterial motility and biofilm formation (51-54). 
Identifying the environmental signals that modulate the 
BswR expression and functionality will answer whether 
BswR represents a previously unappreciated signal 
pathway regulating the virulence and survival capability 
of P. aeruginosa. 
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